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Abstract 

Self-dual string cosmological models provide an effective example of bouncing solu- 
tions where a phase of accelerated contraction smoothly evolves into an epoch of decel- 
erated Friedmann-Robertson-Walker expansion dominated by the dilaton. While the 
00 ' transition to the expanding regime occurs at sub-Planckian curvature scales, the Uni- 

^t^ ■ verse emerging after the bounce is cold, with sharply growing gauge coupling. However, 

\^ . since massless gauge bosons (as well as other massless fields) are super-adiabatically 

^~P , amplified, the energy density of the maximally amplified modes re-entering the horizon 

f~^ I after the bounce can efficiently heat the Universe. As a consequence the gauge coupling 

reaches a constant value, which can still be perturbative. 

> 
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1 Formulation of the problem 

One of the inspiring symmetries of the pre-big bang scenario is the scale factor duahty 
(SFD), ahowing the connection of different cosmological solutions of the low-energy string 
effective action [1]. For instance, in the string frame description, accelerated solutions 
with growing curvature are connected via SFD to decelerated Friedmann-Robertson- Walker 
(FRW) models with decreasing curvature and decreasing dilaton energy density. The same 
solutions can also be described in the so-called Einstein frame, where the dilaton is not di- 
rectly coupled to the Einstein-Hilbert term [2] . In this case a solution exhibiting accelerated 
contraction evolves into a decelerated FRW expansion. 

The connection among different solutions provided by the combined action of SFD and 
time-reversal does not imply an analytical connection of the duality-related solutions, as 
clearly noted in the early days of the pre-big bang cosmology [3]. In spite of this, it is 
known that solutions connecting smoothly two duality-related regimes exist in the pres- 
ence of dilaton potentials which are explicitly invariant under the SFD symmetry [3, 4, 5]. 
Recently its was shown [6, 7] that these models can be derived from a generally covari- 
ant action supplemented by a non-local dilaton potential. In the absence of fluid sources 
and in the Einstein frame metric, these solutions smoothly interpolate between a phase of 
accelerated contraction valid for negative values of the conformal-time coordinate and a 
phase of decelerated expansion valid for positive values of the conformal time coordinate. 
The intermediate high-curvature phase is regular and can be interpreted as cold bounce 
(CB) since the Universe is never dominated by radiation, be it before or after the bounce. 
Furthermore, the nature of CB solutions also implies that the dilaton field never becomes 
constant for ij -^ +oo but is always growing logarithmically. 

A possible way out is to consider the effects of the high-frequency (small-scale) modes 
amplified during the epoch preceding the bounce. In the low-energy effective action also 
Abelian gauge bosons may be present and since their coupling to the dilaton field breaks 
conformal invariance [8] the various modes of the "photon" field can be super-adiabatically 
amplified. While inside the horizon the quantum mechanical fluctuations related to these 
modes are adiabatically damped. On the contrary, as the modes get outside the horizon, 
their energy density is super-adiabatically amplified, thanks to the coupling of the kinetic 
term of the gauge bosons to the dilaton field. When the fluctuations re-enter the horizon, 
after the bounce, their energy density red-shifts as radiation while the energy density of 
the background, still dominated by the dilaton field, red-shifts as stiff matter. Since the 
energy density of the background decreases faster than radiation, there will be a moment 
where the energy density of the radiation background will become dominant, providing a 
mechanism for the gravitational heating of the cold bounce [9]. Well after the bounce the 



dilaton potential is negligible with respect to the dilaton kinetic term; this means that when 
radiation starts dominating the dilaton field will be driven to a constant asymptotic value 
for r] -^ +00. 

The idea that back-reaction effects can produce dynamically a radiation-dominated 
background has been previously discussed within different approaches, but not specifically 
in connection with the problem of dilaton stabilization. The idea of a possible gravita- 
tional reheating of the Universe was pointed out by Parker [10, 11] and also, in a related 
perspective, by Grishchuk [12]. Lukash and Starobinsky [13] considered similar effects but 
paid particular attention to the process of isotropization of anisotropic solutions by particle 
production. In more recent times, Ford [14] discussed the problem of production of slightly 
non-conformally coupled fields in the case of inflationary Universes. Even more recently 
Peebles and Vilenkin [15] used the mechanism of gravitational reheating in the case of 
quintessential inflation (see also [16]). 

In the following the different steps outlined in the previous paragraphs looked into. 
In Section 2 the considerations leading to cold bounces are discussed both in the String 
and in the Einstein frame. In Section 3 the amplification of Abelian gauge bosons are 
accurately computed. Section 4 describes the gravitational reheating of the cold bounces. 
The concluding remarks are presented in Section 5. The Appendices contain useful technical 
results supporting and complementing the results of Sections 2 and 4. 

2 Cold bounces in the string and Einstein frames 

The string frame action leading to cold-bounce solutions [6, 7] can be written as 

S = -^ J <fxsrGe"P[R + G^^d^^dp^ + Vm + S^ (2.1) 

where Xg is the string length scale and the metric signature is mostly minus ^ . In Eq. (2.1) 
Sra account for the possible contribution of matter sources. 

The potential U(^(x)), a local function of ^, is instead a non-local function (yet a scalar 
under general coordinate transformations) of the dilaton owing to the definitions: 

F = U(e-^), e-^(-)= f^^-G{y) e'^^y^ ^G^^-dMy)dMy) ^ {v{^) " viv)) ■ (2-2) 



From the metric G^jy, the induced metric 



dfj^ifdyip 



lilU — G^y (dip)'^ ' 



^ Greek indices run over the four space-time dimensions, while Latin indices take values over the three- 
dimensional spatial geometry. 



and the induced Laplacian 

VV = Tm^V^VV, (2.4) 

can be easily defined. Thus, in terms of Eqs. (2.3) and (2.4), the variation of the action 
(2.1) with respect to G^y and ip leads, respectively, to the following two equations: 



(2.5) 
R + 2G^Pv^Vp^ - G^^da^dp^ + V-^ + e-^^=^h - e'^V'h = 0, (2.6) 

With standard notation, T^ is the energy-momentum tensor of the sources, while the 
integrals 

h = ^Jd^y^-Giy) y'(e-^(^)) 6{^{x) - v{y)) , 

/2 = ^ / A\/-G(y) ^dMy)d^^{y) 6' {^{x) - ^{y)) , (2.7) 

arise directly from the standard variational procedure. 

For a homogeneous, isotropic and spatially flat background, 

500 = 1, g,, = -al{t)5,j, ^ = ^{t), T^ = Ps{t), T^ = -psit)6i, (2.8) 

the evolution equations (2.5) and (2.6) greatly simplify and become perfectly local in time. 
In particular, from Eq. (2.2) 

e"^ = e-^a^, (2.9) 

having absorbed into if the dimensionless constant — ln(/ d'^y/A^), associated with the (fi- 
nite) comoving spatial volume. In the class of backgrounds specified by Eq. (2.8), the time 
and space components of Eq. (2.5), and the dilaton equation (2.6), lead to the following set 
of equations 

^2 _ 3^2 _y^ g^-^^ (2.10) 

Hs = TpHs + ^e^p„ (2.11) 

2ip-Tp^ -m1 + V- — =Q, (2.12) 

■p, + ^Hp, = 0, (2.13) 

where the dot denotes derivation with respect to the cosmic time in the string frame and 
where Hs = cis/as- In Eqs. (2.10)-(2.13) the following quantities have been defined: 

Ps = alps, Ps = alps- (2.14) 



In the string frame (and in the absence of fluid sources) the system of Eqs. (2.10)-(2.12) 
admits various cold-bounce solutions and some examples are provided in the Appendix (see 
also [7]). An interesting case is given, for instance, by the following solution: 

V{^) = -Vbe^^, (2.15) 

_-|l/v^ 



as{t) = ao 



r+ Vr^ + l 
1 



(2.16) 



^ = (/? - 3 In a,(t) = -^ log (1 + r^) + (^q, (2.17) 

where 

t e~^'^o 

The value of t^ sets the maximal value of the Hubble barameter, Hs{t), at the bounce. 
The value of ^pQ sets the typical value of the dilaton at the bounce. 

Using Eq. (2.16), it can be checked that as{t) = aj^{—t) as implied by the invariance 
of the solution under scale factor duality. Other duality-invariant solutions of the system 
of Eqs. (2.10)-(2.13) are reported in Appendix together with a more extensive discussion 
of the cold-bounce solution. 

The conformal time coordinate is the same in both Einstein and string frames [2]. 
Thus, in order to swiftly perform the correct transition from string to Einstein frame, it is 
appropriate to rewrite the system of Eqs. (2.10)-(2.13) in the conformal time coordinate, 
i.e. 

as{r]s)dr]s = dtg. (2.19) 

After this coordinate change, Eqs. (2.10)-(2.13) can be written as: 

ip'^ + en'i - GHsp' = Va^, + a'^e'^ps, (2.20) 

n', = Hsp' - 2nl + ^alps, (2.21) 

dV 
2ip" + ATisif' - en's - 6^2 + y - ^— = 0, (2.22) 

Of 

p', + SHsiPs + Ps) = 0, (2.23) 

where Tig = a'g/^s and the prime denotes derivation with respect to r]s. As already men- 
tioned, the dilaton and the conformal time coordinate do not change in the transition from 
Einstein to string frames 

r]e = ils = r], c^g = c^s = ip, (2.24) 

where r] and ip denote the common values of the conformal time coordinate and of the 
dilaton in both frames. The scale factor, the Hubble parameter, and the energy and pressure 



densities do change in the transition from string to Einstein frame: 

a, = e^/2^, ns = n + ^, p, = e-2V p, = e-2'/'p, (2.25) 

where a, 7i, p and p are the Einstein frame quantities while (p and r/ wih be, respectively, 
the common values of the dilaton and of the conformal time coordinate in both frames. 

Applying the transformations given in (2.25), Eqs. (2.20) and (2.23) lead to the Hamil- 
tonian constraint and to the conservation equation: 

6W2 = pa^ + ^V9'Ve'^Fa2. (2.26) 

p' + m{p + p)-^{p-^p) = Q. (2.27) 

Equations (2.21) and (2.22) lead, respectively, to 

V?" + 2np' + 2n' + ATi^ = pa^, (2.28) 



/2 

if" + 2W(^' + ^ - F + 6H' + 67^2 - ^ = 0. (2.29) 



dV ,. _.,, , ,.,2 V'' 



By linear combination of Eqs. (2.28) and (2.29), the two remaining equations of the system 
can be obtained, namely: 

AH' + 2^2 = _pa2 _ f ^ _ e^PvJ] - e^^a^ (2.30) 

V 2 ) dip 

p" + 2np' + \{p- 'ip)a^ + e'^ {Va" - 1 ^a^") = 0. (2.31) 

Equations (2.26)-(2.27), and (2.30)-(2.31) form a closed set of equations whose solutions 
will now be analysed first in the absence of fluid sources. 

The Einstein frame equations can also be obtained, in generic d spatial dimensions, by 
transforming the string frame action into the Einstein frame and by doing the functional 
variation directly in the Einstein frame variables [7]. Clearly the two approaches lead to 
the same result expressed in eqs. (2.26)-(2.27) and (2.30)-(2.31). 

In the Einstein frame, the fluid evolution and the dilaton are directly coupled in the con- 
servation equation (2.13). In the particular case of radiation, since the energy-momentum 
tensor is traceless, the coupling to the dilaton disappears from Eq. (2.13). 

The solution presented in Eqs. (2.15)-(2.17) cannot be written in simple analytical 
terms in the Einstein frame, but Eqs. (2.26) and (2.27) can be integrated numerically. In 
the Einstein frame description, the asymptotics of the solution (2.15)-(2.17) are 



«<'>-"-f5. «-^-'V^ 



W- = 1^, V'- = -—, (2^32) 



for Tj -^ — oo, and 



a{ri)^a^J^, a+ = e-o/2. / 2(^/3 - 1) 



ip^ =(fO 



2%' ^ V V3 



W+ = T^, f'+ = —, (2.33) 

for ?7 ^ +00. The branch of the solution denoted by minus describes an accelerated 
contraction, since the first derivative of the scale factor is negative while the second is 
positive. The branch of the solution denoted with plus describes, in the Einstein frame, 
a decelerated expansion, since the first derivative of the scale factor is positive while the 
derivative is negative. In both branches the dilaton grows and its derivative is always 
positive-definite (i.e. (p'^ > ) as required by the present approach to bouncing solutions. 

Recalling that, in the Einstein frame, the potential term of Eq. (2.15) and its derivative 
with respect to Ip are 



^e^a^ = We^a^, V = -^0-7^, (2.34) 



riV f.-2'/' 

^e^a2 = We^a\ V = -V^—, 
dip a^ 

Eqs. (2.26)-(2.31) can be numerically integrated across the bounce. 

By setting initial conditions on the asymptotic solution, the integration of the non-linear 
problem can be performed by checking that the constraint of Eq. (2.26) is always accurately 
satisfied. For instance, in the case of the numerical solutions illustrated in Fig. 2.1, the 
constraint is satisfied with a precision of 10~^ over the whole range of the solution. 

In the light of the subsequent applications, it is useful to parametrize the numerical 
evolution of the solution in terms of the gauge coupling and of its derivatives 

g{v) = e^/\ ^'{^) = tL. (2.35) 

9 

In order to find cold-bounce solutions, the system of Eqs. (2.26)-(2.31) is integrated 
in the absence of sources. The initial conditions are chosen in such a way that the gauge 
coupling and the curvature scale, in Planck units, are always minute. Defining r/j as the 
initial time of integration 

., = .ta)«l. ^ = ^«l. (2.36) 



The Haniiltonian constraint of Eq. (2.26) is enforced on the initial data and the initial 
conditions for the scale factor coincide with the solution (2.32). One important parameter 
of the solutions is the width of the bounce, which will be conventionally denoted by tjq. The 
specific value of tjq for a given solution can be determined numerically. Typical values of r/o 
range between rjp and lOOrjp. Furthermore, it turns out that, numerically: 

Hmax ^ —, V'max = —■ (2-37) 

Vo Vo 

These expessions are rather accurately verified for all the sets of initial conditions analysed in 
the present investigation. In the following numerical examples the initial time of integration 
will always be denoted by r]i. For each integration the initial value of the gauge coupling, 
Qi will also be specified. 

In Fig. 2.1 the results of the numerical integration are illustrated, for different sets of 
initial conditions in terms of the evolution of the dilaton (i.e. twice the Neperian logarithm^ 
of the gauge coupling) and of the scale factor. In Fig. 2.2 the evolution of the derivative of 
the dilaton and of the Hubble parameter is reported in conformal time. 

In the left plot of Fig. 2.1, the numerical solution for the evolution of the scale factor 
is compared with the asymptotic solutions provided by the analytical expressions derived 
in Eqs. (2.32) and (2.33). The analytical approximation reproduces rather well with the 
numerical result. The reason for this agreement stems from the right plot in Fig. 2.2 
where the logarithm (in ten basis) of the potential terms is compared with the logarithm 
of the dilaton kinetic energy. Far from the bouncing region the potential is always orders 
of magnitude smaller than the dilaton kinetic energy. Hence the asymptotic solutions are, 
in practice, a good approximation to the full solution not only for rj -^ ibcx) but also, in 
less restrictive terms, away from the bounce. The solutions of Fig. 2.1 and of the following 
plots will be presented in terms of the ratio r]/7]p where rjp = \/^/Mp. In natural Planck 
units, levrG = 2/M| = 1 and rjp = 1. 

The solutions illustrated so far have realistic features. The most notable feature being 
the possibility of describing the transition between the pre-big bang and the post-big bang 
at relatively small curvature. However, various difficulties still have to be faced. The 
first problem is that the Universe emerging after the bounce is cold. The solution, even if 
expanding, is not dominated by radiation as one would like it to be at some later stage of 
the life of the Universe. The second problem is connected with the evolution of the gauge 
coupling. In order to trust our perturbative extimates, it should happen that (7(93) < 1 
throughout all the stages of the model with an aymptotically (constant) value of the order 



^In the present paper the logarithm in ten basis will be denoted by log while the Neperian logarithm will 
be denoted by In. 
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Figure 2.1: In the left plot the result of the integration is illustrated in terms of the (Nepe- 
rian) logarithm of the gauge coupling. In the right plot the scale factor obtained numerically 
(solid line) is compared with the analytical (asymptotic) solutions (dashed lines) following 
from Eqs. (2.32) and (2.33). For both integrations rji = —W^7]p. 





Figure 2.2: In the left plot cold-bounce solutions are illustrated in terms of a' /a and in 
terms of g' /g. In the right plot the dilaton potential is compared with the dilaton kinetic 
term. The various quantities are indicated above or below each curve. As in Fig. 2.1, 
rji = -10"'^7?p. 



of 0.1 or 0.01. In the present case, on the contrary, the dilaton couphng is always growing 
as a power of the conformal time after the bounce. 

3 Amplification of Abelian gauge bosons 



In the low-energy string effective action, massless Abelian gauge bosons are expected [20] 
as a result of the compactification procedure. The modes of the Abelian gauge bosons are 
super-adibatically amplified [8] and they re-enter, at different conformal times, after the 
bounce. The modes of comoving frequencies comparable with the typical curvature at the 
bounce will be the first ones to re-enter. Furthermore the high-frequency modes are also 
the most energetic and effectively behave like a fluid of high-frequency primordial photons. 

Right at the bounce the energy density of the primordial photons is still small. However, 
as the geometry evolves, the energy density of the radiation fluid decreases slower than the 
energy density of the dilatonic sources and, therefore, after a transient time, a radiation- 
dominated phase sets in. The typical transient time depends upon the maximal-curvature 
scale at the bounce and upon the number of massless species, which are parametrically 
amplified. Since in the present case only one Abelian gauge boson is present, the obtained 
estimates represent the minimal achievable efficiency of the process. 

More specifically, consider the case where the coupling of the four-dimensional dilaton 
field to the gauge kinetic term is parametrized as 



1 



3.1) 



Defining the canonically normalized vector potential Ai = g{rj)Ai, in the radiation gauge 



A = 0, 



V-^ = 0, 



3.2) 



the canonical Lagrangian density can be obtained from Eq. ( 3.1) by dropping total deriva- 
tives: 



^{x, V) = ^ 



i^r,A^r + ig-r9^i-i^^A, 



3.3) 



The two physical polarizations of the photon can be quantized according to the standard 
rules of quantization in the radiation gauge in curved space-times [17]: 



and 



A(x,r?) =Y^ 



Mx,v) = Y1 



(27r)3/2 

Sk 

(27r)3/2 



ik-x ^ "t 



dk,c.e'^Ak{vy + flfe ^(^tM'nTe 



\* ik-x 



ik-x I "t „ai 



aMernfc(r?)e-*^ + a^ ^e?nfe(r,)V 



\-k ik-x 



3.4) 



3.5) 



where ef{k) describe the polarizations of the photon and 

The evolution equation for the mode functions will then be, in Fourier space, 



3.6) 



-4'' + 



-i\// 



k'-9{9~') 



A = 0, 



3.7) 



where, following the terminology of [18, 19] the "pump field" ^ can also be expressed as: 

./2 



gig 



-i\// 



V 



2 



(3.8) 



Equation ( 3.7) tells us that all the modes k^ < \g{g ^)"\ are super-adiabatically amplified. 
Naively, the maximal amplified frequency will then be fc^ax ~ \g{g^^)"\- Since, for rj -^ ±00, 



gig 



-iv 



~ 7] 



-2 



modes with krj > 1 are said to be inside the horizon, while modes krj < 1 



are outside the horizon. This not fully accurate terminology will also be employed in the 
following only for sake of simplicity. The maximally amplified modes are then the ones for 
which 



^max — \g[g 



-In// I 



3.9) 



as is illustrated directly in Fig. 3.3 for a particular choice of the parameters. 

The Fourier modes appearing in Eq. ( 3.7) have to be normalized while they are inside 
the horizon for large and negative ry. In this limit the initial conditions provided by quantum 
mechanics are 

1 



Akiv) 
^kiv) 



-ikri 



2k 



-i\ —e 
V 2 



-ikrj 



i 3.10) 



In the limit r] -^ +CX) the positive and negative frequency modes will be mixed, so that 
the solution will be represented in the plane wave orthonormal basis as 






1 



c+ik)e-'^'^ + c.ik)e 



ikrj 



-V2 



c+(A;)e~ 



ikr] 



c-ik)e 



ikrj 



i 3.11) 



^Notice that the "pump field", i.e. the field transferring energy from the background to the fiuctuations, 
is determined by the first and second derivatives of the gauge coupling. 
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Figure 3.3: The pump field for the primordial photons (vertical axis) is plotted. For different 
values of r/o the overall dimensions of the plot are rescaled (since g{g^^)" '^ rjQ ), but the 
relation determining k^axVo always stays the same. In all the plots of the present section 
Vi = -W^VP- 

where c±(k) are the (constant) mixing coefficients. The following two relations fully de- 
termine the square modulus of each of the two mixing coefficients in terms of the complex 
wave- functions obeying Eq. ( 3.7): 



m' 
m' 



\c^{k)\^ = i{AtUk 



AkU* 



|c-(fc)|^ = ^ 



2| A |2 



Uk\' + k'\Ak\ 



( 3.12) 
( 3.13) 



After having numerically computed the time evolution of the properly normalized mode 
functions, Eqs. ( 3.12) and ( 3.13) can be used to infer the value of the relevant mixing 
coefficient (i.e. C-{k)). Equation ( 3.12) is, in fact, the Wronskian of the solutions. If the 
second-order differential equation is written in the form ( 3.7), the Wronskian is always 
conserved throughout the time evolution of the system. Since, from Eq. ( 3.10), the 
Wronskian is equal to 1 initially, it will be equal to 1 all along the time evolution. Thus, 
from Eq. ( 3.12) |c+(A;)p = |c_(A;)p -|- 1. The fact that the Wronskian must always be equal 
to 1 is the measure of the precision of the algorithm. 

In Fig. 3.4 the numerical calculation of the spectrum is illustrated for different values 
of krjQ. In the left plot the mixing coefficients are reported for modes k <^ A;max- In the 
right plot the mixing coefficients are reported for modes around Ajmax- Clearly, from the 
left plot a smaller k leads to a larger mixing coefficient which means that the spectrum is 



11 



rather blue. Furthermore by comparing the amphfication of different modes it is easy to 
infer that the scahng law is |c-|-(A;)p + |c_(A;)p oc (A:/A;inax) "® ) with Ug ~ 3.46, which is in 
excellent agreement with the analytical determination of the mixing coefficients leading to 
Ug = 2^/3 ~ 3.46 [see below, Eq. ( 3.23)]. 

The second piece information that can be drawn from Fig. 3.4 concerns femax; whose 

specific value 

_ \/5-0.5 

can be determined numerically for different values of r/o and will be important for an accurate 
estimate of the back-reaction effects. 

A reasonably accurate estimate of the maximal amplified frequency is necessary in or- 
der to match the analytical estimates with the numerical calculations of the back-reaction 
effects. For the value of fcmax reported in Eq. ( 3.14), the obtained mixing coefficient is 1, 
i.e. |c_(A;inax)| — 1- According to Fig. 3.4 (right plot) as we move from /cmax to larger fc, 
(|c_|_(fc)p + |c_(A;)p) ~ (|c+(/c)p — |c_(fc)p) implying that |c_(A;)| ~ 0. Moreover, from the 
left plot of Fig. 3.4 it can be appreciated that 

|c-(A:^ax)|' = l, log(|c+(A:,nax)|' + |c-(femax)|') = log3~ 0.477. (3.15) 

Thus the absolute normalization and slope of the relevant mixing coefficient can be numer- 
ically determined to be 

k ^-2v^ 



|c-(fc)|^=T— . (3.16) 

\ "'max / 

In Fig. 3.5 (left plot), the predictions of Eq. ( 3.16) at different k modes are compared with 
the numerical calculation. It can be concluded that Eq. ( 3.16) is rather accurate as far as 
both the slope and the absolute normalization are concerned. 

The numerical estimates presented so far can be also corroborated by the usual analyt- 
ical treatment based on the matching of the solutions for the mode functions before and 
after the bounce. The evolution of the modes described by Eq. ( 3.7) can be approximately 
determined from the exact asymptotic solutions given in Eqs. (2.32) and (2.33), and im- 
plying that 99j_ ~ ibv3/ry. Thus the solutions of Eq. ( 3.7) can be obtained in the two 
asymptotic regimes 



AuA^) = ^e^^(-^^y'^H(^\-kr,), r? < -r,,, 



2 



r{i.2) 



_F«(A;r/) + c+e-^"('^+^/2)^(2)(^^) 



, 7] > -r/i( 3.17) 



where Ha ' are Hankel functions of first and second kind [21] whose related indices are 

i^ = — - — , n= — - — . (3.18) 
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Figure 3.4: In the left plot the amplification of modes k"^ <^ |5'(5'^^)"I is illustrated. The 
time evolution of the pump field g{g~^)" was described in Fig. 3.3. In the right plot the 
largest amplified k, i.e. k ~ k^nax, are discussed. The specific value of /cmax discussed in Eq. 
( 3.14) leads to |c_(A;max) — 1- 



The time scale r/i defines the width of the bounce and, typically, rji ~ r/o- 
Phases appearing in Eq. ( 3.17) are carefully chosen so that 



lim Ak 



-ikr] 



Using then the appropriate matching conditions 

■^k-{-Vi) = A,+(^i), 

and defining xi = krji, the obtained mixing coefficients are 



c+{k) 



i!La;^e^-(-+/^+i)/2 



Xl 



Hl^\x,)Hi'\x^) 



+Hl^\x,)Hi'l{x,)+H^'Ux,)Hl'\x,) 



■K 



c_ = i—xie 



'',piA^-t^)/-i 



V + jl + l 

Xi 



H^;^\xi)H^^\x,) 



-H';^\x,)H^^l{x,) + H%{x,)H';^\xr) 



( 3.19) 



( 3.20) 



( 3.21) 



( 3.22) 



satisfying the exact Wronskian normalization condition |c+(A;)p — |c_(A;)p = 1. In the small 
argument limit, i.e. kr]i ~ kr]Q <^ 1 the leading term in Eq. ( 3.22) leads to 



c-{k) 



47r 



,i7r(z.-;.)/2^^^-^(^ + ^ - l)T{fi)T{u) 



( 3.23) 
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Figure 3.5: In the left plot the mixing coefficient computed from Eq. ( 3.16) (crosses) are 
compared with the numerical results (stars). The difference is barely visible. In the right 
plot the logarithmic energy spectrum is computed for k = k^ax and for different sets of 
parameters. The quantity ,^0 = T~(-o/Mp has been introduced. 

If we now insert the values given in Eq. ( 3.18) it turns out that c^{k) ~ 0.41 \kr]i\^^^. 
The spectral slope agrees with the numerical estimate, as already stressed. The absolute 
normalization cannot be determined from Eq. ( 3.23), where the small argument limit has 
already been taken. In order to determine the absolute normalization the specific value of 
kmaxi]! iia,s to be inserted in Eq. ( 3.22). The result of this procedure, taking r]i ~ 770 is 
|c_(fcinax)P = 0.14, which is roughly a factor of 10 smaller than the interpolating formula 
given in Eq. ( 3.16). 

4 Gravitational reheating of a cold bounce 

A reasonably accurate determination of the energy density of the modes re-entering the 
horizon after the bounce and modifying the dynamics of the background can now be pre- 
sented. The expression of the canonical energy-momentum tensor derived from the action 
( 3.1) is 



rpu 



_7? 7?*^" 



4 "/5 



(4.1) 



From the explicit expression of Eq. (4.1) in terms of the gauge potentials Ai, it is possible to 
derive the expectation values of the various components of the canonical energy-momentum 
tensor. While the derivation of these results is included in the Appendix, the relevant results. 
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for the present calculation, are 



1 



Pr = (To°(r?)) = ^ dlnkpik,rj) 






p{k,v) = j^\\^k\^ + 



9 



9 



UlAk + ^kAl 



k-^k 



and 



nJf^W _ _^^i^ 



(r/(r/)) 



{T?iri))=0. 



(4.2) 
(4.3) 

(4.4) 



Sometimes p{k, rj) is called the logarithmic energy spectrum, since it gives the energy density 
of the field per logarithmic interval of frequency. 

For r] -^ — cx) and \krj\ > 1 (i.e. when the modes are inside the horizon at the onset 
of the dynamical evolution) the exact mode functions given in Eq. ( 3.17) imply that the 
logarithmic energy spectrum is 



/j4 



2 + 



^2^2 



+ 



1 



(4.5) 



The first term is nothing but the energy density of the initial vacuum state that can be 
subtracted. The remaining terms are suppressed by powers of Ifcr/I^-*^ and are negligible in 
this regime. 

For \krj\ < 1 the modes are super-adiabatically amplified. In this case the logarithmic 
energy density becomes, from Eqs. ( 3.17) and (4.3) 



pik^v) = ^(-kv) 



rW 



HiXi{-kvr + m'>{-kv)\ 



w, 



(4.6) 



Expanding for \kr]\ < 1, and recalling that Eq. (4.6) is valid for negative r], the leading term 
is 

p{k,v) = ^T\u + l)2''^+\-kvr"'-\ (4.7) 

where 2i/ + 1 = \/3. 



Finally, after the bounce and when the modes have re-entered the horizon, the logarith- 
mic energy density becomes 



P(.k,v) = ^ 



+ 



k^f] 



2|c+(fc)|2 + 2|c_(A:)|2 
(|c+(A;)|2 + |c_(A;)|2 +c;(fc)c„(A:)e2^'='' + clc+e-2i'='7) 



+ii^{cUk)c.{k)e'^''''' - c*_{k)c+{k)e-^'^'^) 
kr] 



(4.8) 



This expression also contains subleading terms going as higher powers of |A;r/| ^. However, 
these terms are negligible when the modes have already re-entered. Hence, using the relation 
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between the mixing coefficients, the leading term appearing in Eq. (4.8) can be rewritten 
as 

p{k,v)c,^\c.{k)\\ (4.9) 

where the initial zero-point energy has already been subtracted. Notice that this final 
estimate is consistent with the interpretation of |c_(A;)p as the mean number of produced 
particles per phase-space interval. 

With this analytical understanding, we are then ready to discuss in detail the effect 
of the produced small-scale fluctuations on the dynamics of the background geometry. As 
expected from the vanishing trace of Eq. (4.1), the effective relation between the components 
of the averaged canonical energy-momentum tensor leads to 3pr = Pr (see also Appendix 
B). Thus the effect of pr as given in Eq. (4.2) should be inserted back into Eqs. (2.26) and 
(2.27) with the result that, for instance, the following equation can be obtained: 

In the numerical code, Eq. (4.10) is consistently solved together with the other equations 
for the dilaton and for the fluctuations of the gauge fields. The problem is simplified by 
the following observation. The integration over k appearing in Eq. (4.2) runs up to fcmax- 
Since the logarithmic energy spectrum decreases at large-distance scales, the integral 

is well estimated by the value of the logarithmic energy spectrum at k^nax- 

pM = ^freo, (4.12) 

where eo = 0.19. 

The result of eq. (4.12) reproduces extremely well the numerical results reported, for 
instance, in Fig. 3.5 (right plot) if we recall the explicit form of fcmax given in the previous 
section. 

The numerical integration may still be rather stiff because of the different scales present 
in the model (in particular the width of the bounce may be much smaller, for instance, than 
the typical time at which radiation becomes dominant). For this purpose the Rosenbrok 
method (particularly suitable for stiff problems [22]) has been used after the bounce. 

Before an illustration of the numerical results, the main semi-analytical expectations will 
be discussed. The typical time scale defining the dominance of radiation can be predicted by 
estimating, numerically, the time of re-entry of k^^ax a-^d by comparing the relative balance 
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between the energy density of the background and the energy density of the fluctuations. 
The time of re-entry of femax is 



11 ^ ^ 

Vmax ~ T ~ —VO, (4-13) 

'I'lnax Cl 



where ei ~ 1.73. 



From this expression, assuming that after the bounce the background is still dominated 
by the kinetic energy of the dilaton, the time rj^, can be estimated to be 



where 62 ~ 12.9. For instance, for Hq/Mp ~ 10~^, and given the indetermination on the 
numerical parameters, ry* ~ 6 x lO^r/p. 

After 77* the evolution is expected to be dominated by radiation, which means that the 
scale factor should go linearly with the conformal time. During radiation the potential 
term becomes even more negligible than in the case of a cold-bounce solution, where it was 
already negligible (recall Fig. 2.2, right plot). Hence, the evolution of the dilaton will be, 
from the background equation (2.31), 

(aV)' = (4.15) 

which also implies that 

ip^ipi-ip2{ — 1, (4.16) 

where cpi and 992 are integration constants. 

A useful function, which can be used as a diagnostic of a fully achieved transition to 
radiation, is 

zin) = ^. (4.17) 

It should be noticed that 2(77) ~ a(rj), if evaluated on the cold-bounce solution, while z{ri) 
is strictly constant , in the case of a radiation-dominated background, by virtue of Eqs. 
(4.15) and (4.16). 

The analytical expectations suggested by Eqs. (4.14)-(4.17) are confirmed by the numer- 
ical integration of the full system, which is illustrated in Figs. 4.6 and 4.7. The constancy 
of z{ri) is exactly what is numerically observed. In Fig. 4.6 (left plot), after a typical time- 
scale 7]^: ~ lO^r/p, z(ri) becomes constant while it should still increase as 2(77) ~ a{ri) ~ y^ 
if the Universe would not be dominated by radiation. At the same time (Fig. 4.6, right 
plot) the scale factor grows linearly and (Fig. 4.7, left plot) the dilaton goes to constant, 
as implied by Eq. (4.16). 
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Figure 4.6: In the left plot the evolution of z{r\)^ defined in Eq. (4.17), is illustrated. This 
simple function is a diagnostic of a fully achieved transition to radiation where it is constant, 
while it should grow as y^ (dotted line) in the case when the back-reaction of photons is 
not included. In the right plot the evolution of the scale factor is illustrated with (solid 
line) and without (dashed line) back-reaction effects. 
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Figure 4.7: In the left plot the evolution of the logarithm of the gauge coupling is illustrated. 
In the right plot the evolution of a' /a is presented. 
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5 Concluding discussions 

Cold-bounce solutions have virtues and problems. Their main virtues are that they lead to 
computable models of the bouncing dynamics at low curvatures. The main problems are 
related to the fact that the Universe emerging after the bounce is rather cold and dominated 
by the dilaton. This occurrence also implies that the dilaton will always be growing, before 
and after the bounce. 

In the present paper a model of gravitational heating of the cold bounce has been 
proposed. The model considers only the effect of one species of Abelian gauge bosons, 
which are super-adiabatically amplified. The frequencies that are maximally amplified are 
comparable with the typical curvature of the Universe at the bounce and effectively behave 
like a fluid of primordial photons. The following results have been illustrated: 

• an accurate numerical method for the calculation of the amplification of the primordial 
photons has been developed; 

• taking into account the back-reaction of the primordial photons, the cold-bounce 
solution can be consistently heated; 

• the transition to radiation can occurr for sub-Planckian curvature scales; 

• solutions have been presented where the dilaton goes to a constant value and the 
asymptotic value of the gauge coupling is always smaller than 1. 
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Tp^ -dH"^ -V = e'^p, 






H = ipH + ]^e^p, 






2ip-Tp'^ -dH^ + V - 


dV 

dip 


= 0, 


"p + dHp = 0. 







A Solutions with and without fluid sources 

Equations (2.10)-(2.13) can be written in generic d spatial dimensions as: 

(A.l) 
(A.2) 

(A.3) 
(A.4) 

In this Appendix, cold-bounce solutions with and without fluid sources will be presented. 
It will be shown that the naive addition of fluid sources does not stabilize the evolution of 
the dilaton field after the bounce. Since throughout this Appendix we will always consider 
solutions in the string frame, the subscript has been omitted in the quantities appearing in 
Eqs. (A.1)-(A.4). 

A.l Solutions without fluid sources 

When fluid sources are absent there is a class of regular solutions in the presence of ex- 
ponential potential of the form V^(^) = — Voe"*^. Defining the appropriate variables of the 

phase space 

Tp ^f-V 

u = —= — , V = —= — , (A. 5) 

VdH VdH 

the constraint (A.l) becomes u'^ + v'^ = 1 and the equations for u and v are 

^ = V3(l - u') - ^■ylv' s -^(a - 2)(1 - „=), (A.6) 

where, in the first equation, the equality follows from u'^ + v'^ = 1 and where r = log a. Eqs. 
(A.6) and (A. 7) then are in the form of an autonomous system in the plane {u,v). From 
Eq. (A. 5) it is clear that 

t -^ —t, u ^ u, V ^ —V. (A. 8) 

Hence, the system (A.6)-(A.7) is symmetric for u ^ u and v -^ —v when r -^ — r. Thanks 
to the above symmetry the phase space of the autonomous system (A.6)-(A.7) is fully 
characterized by half of the unit disk, for instance, the one for v > 0. 

The critical points of the system are given by 

Uc = ±1, Vc = 0. (A.9) 
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The plane autonomous system can then be hnearized around the critical points given in Eq. 
(A. 9), namely 

u = Uc + eu, v = Vc + ey, (A.IO) 

— ^ = 2iJ,Uceu, (A.ll) 

(IT 

-p = fiUcCv, (A.12) 

where 2fi = \/d{a — 2). If Uc = 1 the critical point is an unstable node since the two 
eigenvalues have the same (positive) sign. If Uc = —1 the critical point is a stable node since 
the two eigenvalues have the same sign and are both negative. The unstable node (1,0) 
corresponds to the accelerated branch of the solution where the dilaton coupling increases 
and the curvature scale also increases. The stable node (—1,0) corresponds to the dual 
(decelerated) solution with decreasing dilaton and decreasing curvature. 

In this parametrization the evolution equations can be easily integrated. From Eq. 
(A. 6), integrating once and inverting the obtained expression: 

<r) = j^^^, (A.13) 

where the integration constant C has to be selected such that the appropriate boundary 
conditions are satisfied. The constant C can be chosen in such a way that u{t) is everywhere 
regular in its domain of definition and, in particular, for r ^ 0. With this in mind, C = — 1 
in order to prevent a possible singularity in r = 0. From Eqs. (A. 7) and (A.13), the full 
result is 

u{t) = -tanh^r, (A. 14) 

v{t) = -^—, (A.15) 

cosh HT 

where the possible further integration constant is fixed by the constraint u^ + v'^ = 1. 

The solution given by Eqs. (A. 14) and (A.15) interpolates between the unstable node 
(1,0) and the stable node (—1,0), as is apparent from the asymptotic behaviour of {u,v) 
for r -^ iboo. Recalling now the explicit definition of u and v in terms of the original 
background fields, i.e. Eqs. (A. 6) and (A. 7), and recalling the explicit parametrization of 
the exponential potential, the following chain of identities can be obtained: 

u(r) = ^^ = ^-^ = -tanh^r, (A.16) 

VdH VddT 

v{t) = "^ 7- = -^—. (A.17) 

y/dH cosh/ir 
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Integrating once Eq. (A. 16) and substituting in (A. 17): 

2 

Tplr) = ifQ Incosh(^r). (A. 19) 

a — 2 

The parametrization in terms of r = In a is then more useful than the one in terms of t. 

In fact, Eqs. (A. 18) and (A. 19) have simple analytical expressions in terms of r, but not 

necessarily in terms of t. Recalling that, by definition, dr/dt = H{t), the explicit relation 

between r and t can be worked out by integrating once H{t) from Eq. (A. 18). For instance, 

in the case a = A 



ip = LfQ — \n cosh v3r, 
dr e^^o^/T^ 



(A.20) 



dt -^/S cosh "v/Sr 

defining to ^ = e2'^o/i^, 

t/to = sinh \/3r, (A.21) 

i.e. r = Ina = (l/\/3) In [(t/to)^ + VW^oP~+^]- This last solution is the one quoted in 
Eqs. (2.15)-(2.17). Note that if we do not want H{t) blowing up for |r| -^ oo, we have 
to demand a > 2. In the case a = 2, /x = and Eqs. (A.6)-(A.7) imply that u and v are 
constant. 

A. 2 Solutions with fluid sources 

If fluid sources are included, various solutions are possible. Consider first the case p = 0. 
In this case Eq. (2.11) can be immediately integrated and an exact regular solution of the 
remaining equations is 

y = -^06^-^16^^, (A.22) 



with 



H ' *° 



1 t^ 

2^^k 

-^ „-d 



f = ^G- -log [-2+1], 



P = Poa-, (A.23) 



subject to the conditions: 



Po = Vc 



0, 



VI e^'^" = \. (A.24) 
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Consider now the case where the barotropic index 

w = ^ = ^ (A.25) 

P P 

is constant, and the dilaton potential is y = —Voe'^^. In this case, Eqs. (A. 2) and (A. 3) 
can be integrated : 

H = w—{x + xi), 

Tp = - — {x + xo), (A.26) 

where 

dx L / . „x 

and xi, Xo are integration constants ^. 

From Eq. (A.26) it can be checked that cp = Tp + 3H is always positive only in the case 
w = 1/3, corresponding to the case of a radiation fluid. In this case, an appropriate choice 
of the parameters xi and xq allows a simple analytical solution, which can be expressed 
directly in terms of x, i.e. the new time variable. If w = 1/3, then Eq. (A.26) can be 
inserted back into Eq. (A.l), with the result 

f,¥ = ^^ ^A 28^ 

3(x + xo)2-(x + xi)2 + 3yo^2' 

Now we can choose constants to further simplify the general solution. Choose, in particular, 

Xo = L, xi = 3x0, VqL^ = -L^. (A.29) 

Then, it can be checked that the explicit solution is, in this case 

a(x) = Vx2 + L2e3^''^*^°(^/^), 
a' X + 3L 



a x^ + L^ ' 
_, _ 3(x + L) 
'^ ~ x2 + L2 ^ 
dx pqL 1 

Hi ~ ~a(x)' 



(A.30) 



The shifted dilaton is 



3 

Tp = LpQ — ?> arctan (x/L) — - In (x + L ). (A. 31) 



'Since p has dimensions of an inverse length, x and t have the same dimensions. 
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Looking at the constraint, we should also bear in mind that 



e-» = I'PoL'. 



(A.32) 



There is a third class of possible solutions of the system of Eqs. (2.10)-(2.13). This 
class corresponds to the case when the barotropic index is a function of time. In this case, 
Eq. (A. 26) is modified, i.e. 



H = — I dxw{x), 



"P 



-j-{x + xo). 



(A.33) 



A particularly simple example of this class of solutions is provided by the profile [2] : 



1 



w{x) = - 



^ ■ ''x2 + X? 



(A.34) 



By integrating the equations of the background with a method similar to the one introduced 
in the case of constant w solutions, analytical solutions can be found in the case of V = 
—Vqc^"^. Notice that in the case described by Eq. (A.34) the barotropic index interpolates 
between w = —1/d and w = 1/d, i.e. for d = 3, w = —1/3 and w = 1/3. The full solution 
is then: 

2 



(p = ipO 



p'PO 



d 



d- 1 



log 



^1 +1 

Xl 



XlJ 

d+1 
2(d-l) 



d-l 



oVa 



+ 1 



1 

" d-l 



(A.35) 



B Averages of the energy-momentum tensor 



In this Appendix the averages of the canonical energy-momentum tensor given in Eq. (4.1) 
will be computed. In the radiation gauge the components of the energy-momentum tensor 
read 



^°"2^ 



{A'if + ( ^\ A} + 2^A'iA' + {diAjf - d^Ajd^A' 



9 



9 



fB.l) 



/\ 2 



T^ = ^(dr,A,d^A^ + 2(^yAA^ + 2^ 



C/y^ ./H.^ ^T? "^ 



2diAkd^A'' + 2dkAid^A'' + ^ 



Cf) ./T.'i (yf) v^ 



yfj^/^iL/fj 
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2^d^AiA' - ( -^ AiA' + diAkd'A'' - diAud'^A' 



.9_ 
9 



' a4 



9 

d^A'' + ^^M ( diAk - dkAi 



(B.2) 
(B.3) 



Note that in Eqs. (B.1)-(B.3) the spatial indices are flat. Furthermore, in Eqs. (B.2) and 
(B.3) the derivation with respect to ij has been indicated with 9^ (and not, as usual, with 
a prime) only in order to avoid confusion with controvariant (flat space) indices. 

Recalling now the expression of the field operators given in Eqs. ( 3.4) and ( 3.5), the 
various expectation values implied by Eqs. (B.1)-(B.3) can be evaluated separately: 



{Ai{x,ri)Aj{y,r])) 

{Ui{x,r])Uj{y,ri)) 
Sk 



(27r)3 

d^k 

(2^ 



P^,{k)\Ak{r]re 



2 -ik-{x-y) 



P^jik)\Ukir])\'e 



2 -ik-{x-y) 



i2nY 



Pi,ik)[Ukir]yAkin) + Ukiv)Akir]y]e 



*]p~ik-{x~y} 



where 



Kirvj 



Using Eqs. (B.4) (for spatially coincident points) into Eqs. (B.l) we obtain 



(T.»(,)> = i/ 



d^k 

(2^ 



\Ilk\' + 



9 



9 



TllAk + TlkAl 



(B.4) 
(B.5) 

(B.6) 



leading, after angular integration, to Eq. (4.2). The same type of calculation can be 
performed in the case of the expectation value of Eq. (B.2). The result in this case can be 
written as 

1 /• d^k 



{Tim 



a* J (27r)3l^'^^^ 



\Iik? + ( -)Vfc|' + - ( AkK + Allik 



-AhknAk? + 5l 



k'-\'-)\Ak\'' 



mk\''--[KAk + ^kAk 



(B.7) 



If f{k, rf) is a generic scalar function of the modulus of the momentum and of the time the 
following identities hold: 



d^kkik^fik,7]) = —S{ I k^f{k,r])dlnk, 
d^kk'f{k,r]) = 0. 



(B.8) 
(B.9) 



Using Eq. (B.8) into Eq. (B.7) and recalling the explicit expression of the projector P/(fc), 
we can easily obtain 

(Ti) = -lp,6i^-p,6i, (B.IO) 
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as can be argued from the fact that the classical energy-momentum tensor is traceless. 
Finally, the expectation value (T^) always vanishes by virtue of Eq. (B.9). 
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